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RESEARCH MEMORANDUM 

DETERMINATION OP MINIMUM MOMENTS OP INERTIA OP 
ARBITRARILY SHARED AREAS, SUCH AS 
HOLLOW TURBINE BLADES 
By Sel Gendler and Donald P. Johnson 


SUMMARY 

A simple accurate method is presented for approximately deter- 
mining the minimum moment of inertia of an arbitrarily shaped area, 
such as the seotion of a hollow turbine blade. The praotieal appli- 
cation of this method involves a simple routine tabular procedure. 
Incidental to finding the minimum moment of inertia, this tabular 
procedure also gives the area, the position of the center of gravity, 
the moment of inertia about any desired axis, the product of inertia, 
and the principal axes of inertia. Two examples are worked out in 
detail: an ellipse tilted at an angle of 30 , and an airfoil seotion 

In general, comparison with the known values for the ellipse 
showed agreement to within 0.5 percent. This method was also used to 
calculate the variation in minimum moment of inertia along the length 
of a typical hollow turbine blade. 


INTRODUCTION 

In several types of problems involving vibrations of an elastic 
body, it is necessary to know the area, the plane product and the 
moment of inertia about a given axis, the minimum moment of inertia, 
and the variations of these quantities along an axis through the body 
The most obvious method, that of "counting squares", can be applied 
to obtain the area, the product of inertia, and the moment of inertia 
about a given axis. The minimum moment of inertia can be approxi- 
mated by a trial -and -err or procedure of passing several axes through 
the center of gravity and calculating the moment of inertia about 
each of these axes. This method is long and laborious and the 
accuracy attained depends a great deal upon the skill of the operator 
A more expedient method, which is in general use, is that of 
numerical integration. 
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As a result of an Investigation, conducted at the NACA Levis 
laboratory, of the vibration problems of compressor and turbine 
blades in general and hollow turbine blades in particular, a simple 
accurate method of calculating the moment of inertia about a given 
axis and the minimum moment of inertia of arbitrarily shaped sec- 
tions has been developed. This method combines the well-known 
theories of moments of inertia and numerical integration to provide 
a simple straightforward tabular procedure for calculating the 
section properties. The area, the position of the center of 
gravity, the moment of inertia about any given axis, and the pro- 
duct of inertia are obtained incidental to finding the minimum 
moment of inertia. The values obtained are dependable and may be 
calculated to any reasonable degree of accuracy without unduly 
Increasing the amount of labor. The procedure is illustrated by 
two simple examples, an ellipse tilted at an angle and an airfoil 
section, and is finally applied to finding the variation of the 
minimum moment of inertia along the length of a twisted and tapered 
hollow turbine blade. 


SYMBOLS 


A 

\ 

C 

c 

F 


I m,r 

*x 


L X.O 


The following symbols are used in this report: 
cross-sectional area 
cross -sectional area of root section 
boundary 
ohord 

product of inertia about the x,y-axes 

product of inertia with respect to axes through centroid 
parallel to x,y-axes 

minimum moment of inertia 

minimum moment of inertia at root 

moment of inertia about x-axis 

moment of inertia about axis through centroid parallel to 
x-axis 

moment of inertia about y-axis 
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I y 0 moment of inertia about aria through centroid parallel to 
’ y-axis 

minimum x inertia ratio 
I m /ly^o minimum y inertia ratio 
Z length of blade 


n 


*,y 

xSy' 


*>7 

*0>7q 

a 

e 


number of stations into -which area is divided 
reference coordinate axes 

principal axes about one of which moment of inertia of area 
is minimum 

location of centroid of area referred to x,y-axes 
location of origin of x^y'-axes referred to x,y-axes 
inertia parameter, 2 F 0 /(l y ^ o - I^ 0 ) 

angle of x* -axis to x-axis measured in positive counter- 
clockwise direction 


METHOD 


Let x,y be any set of reference axes with respect to which 
the coordinates of the area under consideration are known (fig. l). 
Let x’,y' be the principal axes at the point (xq,Jq) , such that 

the moment of inertia of the area is a minimum about the x'-axis, 
and let x,y be the coordinates of the centroid of the area with 
respect to the x,y axes. Then (reference l) : 


Xq = x 
7q - J 



( 1 ) 


2 TV 


I -I 
7,o x,o 


tan 26 


( 2 ) 
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Hence, It can be derived that 


where 


I m - -F 0 cot 6 + Iy^ 0 



o _I x, 


) 




(3) 


a 


2 F. 




(4) 


Equation (l) states the well-known fact that the moment of 
inertia will be least about an axis passing through the centroid. 
Equation (2) gives the angle that this axis makes with the x-axis. 
In equation (3), which gives the minimum moment of inertia, the sign 
before the radical is so chosen that the first term on the right 
side of the equation is negative. 

If the equations of the contour are known, equations (l) to 
(3) can be solved exactly for the desired quantities. (See the 
appendix.) These equations, however, are generally unknown and 
approximate methods must be used. 


Assume that the values of y defining the area boundary are 
known for n stations along the x-axis. These values may be obtained 
by direct measurement from a drawing, or as is generally the case for 
airfoil sections, they may be directly obtained from a table that 
gives the values of y in percentage of chord for a number of stations 
along the chord. The following equations then approximately hold: 
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Also 



1=1 


F * ~ \ **1 > 
1=1 ‘ 

lx - - | j>V 
1=1 




I y = 


-Z 


i=i 


yixi Axi 



(5) 


x x,o - x* - y 2a 

■VfO “ -V “ x A 

F 0 = F - x y A 


The values x^ and y^ should he taken at the midpoint of the 
interval Axj. It is convenient to take equal intervals Ax, because 
by so doing, Ax^ moves to a position in front of the summation sign 
and thus considerably shortens computation. 
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The problem then reduces to one of tabulating n values of x 
and y and obtaining the sums of the various products Indicated In 
equations (5) . The minimum moment of Inertia Is then found by use 
of equations (3) and (4). The accuracy with vhlch the factors In 
equations (5) sure obtained increases vith n, the number of stations 
used. In general. It vlll be seen from the examples that good 
accuracy can be obtained by use of relatively few stations. 


In carrying out the preceding calculations, care must be taken 
as to the sign of the quantity Ax. The summation, as shown by the 
arrows in figures 1 and 2 , is in such a direction that the area 
under consideration is always to the left. If the summation is in 
the direction of Increasing x, the sign of Ax should 
be positive; If the summation Is In the direction of decreasing x, 
the sign of Ax should be negative. Thus, In figure 1, Ax is 
positive for curves and C 3 and negative for curves Cg and 

C 4 . This convention has been followed In the derivation of equa- 
tions (5) and must be followed in their application. 


Graphical solution of 3^. - Equation (3) Is made nondlmenslonal 


If both members are divided by I^g. 



( 6 ) 


Equation ( 6 ) represents an infinite family of straight lines 
with the minimum inertia ratio, Im/li /0 > as ordinate and 
(l y 0 /l x#0 ) - 1 as abscissa. (In actual practice. It will be 
found more convenient to use as bhe abscissa.) All the 

lines have a common intercept at = 1 and the slopes of 

the lines depend on the inertia parameter a as follows: 


Slope = ^ 

It is evident that the abscissa cannot become less than 0 
and the ordinate cannot be less than 0 or greater than 1. More- 
over, there will be no lines with slopes between 0 and 1. Such a 
family is plotted in figure 3. The lines giving the maximum value 


1 + (l+ar) 
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of I x ,o t>© merely the reflections of the lines in figure 3 

about the point (1,1). If the values of a, I^q, and I y Q have 

been calculated for a section, the values of I m can be obtained 
from figure 3. 


E xamp le 1« - Figure 2(a) shows an ellipse tilted at an angle of 

,2 2 

30° to the x-axis. The equation of this ellipse is + y* = 1. 

The area, the location of the centroid, the moment of inertia about 
the x-axis, the principal axes through the centroid, and the minimum 
moment of inertia are to be determined. The ellipse was divided 
into ten intervals. The quantities in equations (5) are calculated as 
shown in table I. The following values were obtained: 


A = 9.516 
x = 0.000 
y » 0.000 


I x = 7.212 
3^ = 2.363 
6 = 0.521 radian 


The exact values for the ellipse are: 


A = 9.425 
x = 0.000 
y - 0.000 


I x = 7.069 
= 2.356 

9 = 0.524 radian 


The decrease in percentage error of 3^. due to choosing a 

larger number of intervals over the ellipse is shown in the follow- 
ing table: 


Error in calculated 
value of Ij^, percent 

10 

2 

.3 

.1 


Number of stations 


2 

4 

10 

20 


•ttTfmvpi <=) 2. - A Clark Y airfoil is shown in figure 2(b). The 
shape of this airfoil in terms of the chord is given in table II . 
The q u a nt ities in equations (5) are obtained as shown in table III. 
The final answers are : 



8 


A = 
x = 

y - 

f = 


815 - 7 (#) 
42 - 05 fe ) 



HACA EM E9H10 



G = 0.0358 radian 


DISCUSSIOH OF RESULTS 

Examples 1 and 2 with tables I and III show an application of 
the method for obtaining the minimum moment of inertia of arbitrarily 
shaped areas. As has already been, pointed out, equations (5) are 
approximate, the degree of approximation depending upon the number 
of stations considered. It is evident that in the oase of the 
ellipse, the error involved in the use of ten stations Is less than 
0.5 percent. Thus, it is unnecessary to consider many stations to 
achieve good accuraoy. 

The number of stations required to give a certain degree of 
accuracy will, of course, depend on the shape of the area under con- 
sideration. In order to obtain good over-all accuracy, areas that 
have sharp peaks in their contours, such as a cross-shaped section, 
will require more stations in the vicinity of the peaks than in the 
more uniform parts. Fortunately, most objects dealt with in physical 
problems do not have such sharp peaks. If such an object is being 
considered, however, care must be taken in selecting the number of 
intervals and in distributing the intervals. It may be advantageous 
to take small intervals for certain parts of the contour and large 
intervals for the rest of the contour. 

Possible sources of error should be pointed out. If the angle 6 

is very close to 45°, the denominator of equation (2) becomes very 

small. If small errors do occur in the calculation of I_ _ and 

y 

I x 0 , it is possible that the wrong sign will be obtained for the 
square root of equation (3). If ls a 8111611 quantity, a large 

error in the value of ^ could be introduced. Such an error could 
easily be prevented because by inspection, it should be obvious 
whether G is 45° or -45°. In any case, it is better to choose the 
x^y'-axes such that G vill not be close to 45°. 
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If the coordinate axes are so chosen that the angle 6 is known 
to he small, an approximate form of equation (3) can he used. For 
small 6, a is small and a? becomes negligible compared to 1. 
Equation (3) then becomes 

^ = I t^o-^oX 1 ± « + **,0 

Thus, 3^ as a minimum will equal either 1^ Q or 3^. Q , 
whichever is smaller. If, however, the absolute value of Q 

is large compared to I x Q , the neglect of even a small quantity 

g } 

a can introduce large errors in 3^. 

The method described was used to calculate the minimum moment of 
inertia of various sections along the length of a hollow turbine 
blade. Five sections, as shown in figure 4, including the tip and 
root section were used. The results and the area variation are 
plotted in dimensionless coordinates, as shown In figure 5. These 
variations must be known in order to calculate accurately the natural 
vibration frequencies of the turbine blades. 


SUMMARY OF RESULTS 

A simple tabular method for obtaining the minimum moment of 
inertia of an arbitrarily shaped section was developed. The area, 
the position of the centroid, the moment of inertia about an arbi- 
trary set of axes, and the product of Inertia were obtained as part 
of this tabular procedure. Although the accuracy obtained using a 
given number of stations depends upon the irregularity of the area 
under consideration, it was shown that, in general, high accuracy 
could be obtained with the use of comparatively few stations. In 
the case of an ellipse tilted at 30?, accuracy of greater than 0.5 
percent was obtained by the use of ten stations. A final applica- 
tion was made in calculating the variation of the minimum moment 
of inertia along the length of a hollow turbine blade. 


Lewis Flight Propulsion Laboratory, 

Rational Advisory Committee for Aeronautics, 
Cleveland, Ohio. 



10 


MCA a* E9H10 


APPENDIX 

EXACT SOLUTION OF MINIMUM MOMENT OF INERTIA 


Minimum moment of Inertia. - Let x,y be the coordinates of a 
point of a closed contour relative to a given set of x,y-axes; let 
x',y ' be the coordinates of this point relative to the principal 
axes of the contour (about which the moment of inertia is a mini- 
noun); and let 9 be the angle between the two sets of axes. Then, 

3^ = cot 9 (xyA-pj + I y -x ^ A (Al) 


where 


tan 29 


2(F-xyA) 


2 F, 


_ 2 


(Iy-x A) - (Dj.-y A) 


C y,o _I x,o 


(A2) 


By taking tan 29 = fL_tan0_ 3^3 letting a = •= ^ f the 

1-tan 9 ■ L y,o“ x i,o 

following equation is obtained by substituting in equation (Al) and 
carrying through the proper trigonometric manipulation: 


= 2 ^ I y,o“ I x,o^ 


|_1 + (1-wf) J + 3^ 


(A3) 


The sign preceding the radical is determined by the condition 
that I^. q-Ijj! and Iy #0 -I m must each exceed zero. Thus, if 

(l_ -I ) > 0, a minus sign is used and if (l_ -I ) > 0, a 
y,o x,o ' y,o 

plus sign is used. 

Exact solution. - If P(x,y) and ^ are continuous 

within and on the boundary of region A, then by use of Gauss' 
theorem (reference 2) 



(A4) 


where the integration on the right is performed in a positive sense 
around the boundary C. 



NACA EM E9H10 


11 


i-* 
CD 
f— ' 


From equation (A4) and the fact that the line integral around 
closed boundary E(x) dx equals zero, the following equations 


are evident: 


Also 


- SS 44 • * St ax 

a c 

-X JJ iaA '- a Z 3 ^ 41 


i SS A 7 “ * - £ f c S 41 


F = 


lx « 


X 7 = 


SS ^^--1 f a 
SS 


xy 2 dx 


y 2 dA = - | jT y 3 dx 
'A 3 J C 


x 2 dA 


r x,o • J x - y A 


- ^xVax y 


(A5) 


I y,0 = Z y - x A 


Fq = F - x y A 

If the equation of the boundary of the section under consider- 
ation is known, the quantities in equation (A5) can then be exactly 
determined by carrying out the line integration indicated. The 
integration is in a direction such that the area is always to the 
left. 
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TABLE I - SOLUTION OP EQUATIONS (5) FOR ELLIPSE 


Ax 

X 

y 

y 2 


XT 2 

*2y 

xy ■ 

-0.52915 

2.3812 

-0.6842 

0.4681 






1.9520 

-.1067 

.0114 






1 . 3229 

.3273 

.1071 






.7937 

.6889 

.4746 






.2646 

.9973 

.9946 






-.2646 

1.2592 

1.5856 






-.7937 

1.4745 

2.1742 






-1.3229 

1.6367 

2.6788 






-1.8520 

1.7263 

2.9801 






-2.3812 

1.6727 

2.7979 





Total 


8.9920 

14.2724 

20.4439 

-15.9532 

16.1181 

-11.4320 

0.52915 

2.3812 

-1.6727 

2.7979 






1.8520 

-1.7263 

2.9801 






1.3229 

-1.6367 

2.6788 






.7937 

-1.4745 

2.1742 






.2646 

-1.2592 

1.5856 






-.2646 

-.9973 

.9946 






-.7937 

-.6889 

.4746 






-1.3229 

-.3273 

.1071 






-1.8520 

.1067 

.0114 






-2.3812 

.6842 

.4681 





Total 


-8.9920 

14.2724 

-20.4439 

15.9532 

-16.1181 

-11.4320 



A * - EyAx = 9.5162 

_ £ xy A x 

x - - — ^ - 0.0000 

y » - £y ^ A Ax « 0.0000 

P - - £ *y 2 . 8.4416 

2 

£y® Ax 

lx ■ " 3 ' 

I y = - Sx2y Ax * 17.0578 


I 


x,o 


^y.o 


a 




e 


= I x - y 2 A = 7.2119 
m Iy - x 2 A * 17.0578 

* P - xf A - 8.4416 

2P 0 

■ =1.7147 

^T.o “ ^“x.o 

■^■.o ■ i x,o r . o.i-i 

* 2 [l-(l+ «*) J+Ix,o - 2 

1 ■* gp 1 

“ — tan -1 Q “ 0.5214 radian 

2 I - I 

r»o x x,o 


3628 
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TABLE III - SOLUTION OP EQUATIONS (5) FOR CLARK Y AIRFOIL 


Ax 

X 

7 

y 2 

y 3 

xy 2 

X 2 y 

xy 

H 

O 

10 

9.60 

92.1600 






20 

11.36 

129.0496 






30 

11.70 

136.8900 






40 

11.40 

129.9600 






50 

10.52 

110.6704 






60 

9.15 

83.7225 






70 

7.35 

54.0225 






80 

5.22 

27.2484 






90 

2.80 

7.8400 





Total 


79.10 

771.5634 

7925.4640 

30,051.6090 

185,617.0000 

3389 . 3000 

-5 

2.5 

6.500 

42.2500 






97.5 

.805 

.6480 





Total 


7.305 

42.8980 

275.1466 

168.8050 

7693.1563 

94.7375 

10 

10 

0.42 

.1764 

• 





20 

.03 

.0009 






30 

0 

0 






40 

0 

0 






50 

0 

0 






60 

0 

0 






70 

0 

0 






80 

0 

0 






90 

0 

0 





Total 


.45 

.1773 

0.0741 

1.7820 

54.0000 

4.8000 

5 

2.5 

1.47 

2.1609 






97.5 

0 

0 





Total 


1.47 

2.1609 

3.1765 

5 . 4022 

9.1875 

3.6750 



A = 


/ C 2 X 

I,* - 2y3 £lX - 26, 870 ( 

° 4 ) 

O 

H 

x 3 V 

< 10 s J 


_ Sxy Ax 
” A 


42.05 


£y Ax = 4.853 

2A \ 10 2 


F = - 


Zxy 2 Ax _ 


10 2 J 

te) 

C 4 


150,600 ( — ^ 

\ 10 8 J 


I v = -2x2y Ax - 1,894,000 
* \108/ 


a * 0.07169 


• 7087 (jss) 


0 = 0.0358 radian 
































(b) Clark Y airfoil. 

Figure 2. - Areas under consideration for minimum 
moment of inertia. 




Minimum x inertia ratio 



Minimum y inertia ratio I Xf0 

Figure 3. - Variation of minimum x inertia ratio I]„/l x#0 with ratio of minimum x inertia ratio I,n/ X>0 

to minimum y inertia ratio Iy f0 /*x,o" 
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Figure 4. - Variation of prpss secti on of hol low turbine 
blade. Length, 4 inches. 





